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ABSTRACT 

The surprising discovery by MAGIC of an intense, rapidly varying very high energy {E > 50 
GeV) emission from the flat spectrum radio quasar PKS 1222+216 represents a challenge for 
all interpretative scenarios. Indeed, in order to avoid absorption of 7 rays in the dense ultravi- 
olet radiation field of the broad line region (BLR), one is forced to invoke the existence of a 
very compact (r ^ 10^* cm) emitting region at a large distance {R > 10^^ cm) from the jet 
base. We present a scenario based on the standard blazar model for PKS 1222+216 where 7 
rays are produced close to the central engine, but we add the new assumption that inside the 
source photons can oscillate into axion-like particles, which are a generic prediction of many 
extensions of the Standard Model of elementary particle interactions. As a result, a consider- 
able fraction of photons can escape absorption from the BLR much in the same way as they 
largely avoid absorption from extragalactic background light when propagating over cosmic 
distances. We show that observations can be explained in this way for reasonable values of the 
model parameters, and in particular we find it quite remarkable that the most favourable value 
of photon- ALP coupling happens to be the same in both situations. An independent laboratory 
check of our proposal can be performed by the planned upgrade of the ALPS experiment at 
DESY. 

Key words: radiation mechanisms: non-thermal — 7-rays: theory — galaxies: individual: 
PKS 1222+216 



. 1 INTRODUCTION 

. Blazars dominate the extragalactic 7-ray sky, both at high energy 

■ (> 100 MeV) and at very high energy (VHE) (> 50 GeV). Their 
' powerful non-thermal emission, spanning the entire electromag- 

■ netic spectrum, is produced in a relativistic jet pointing toward the 

■ Earth. The spectral energy distribution (SED) shows two well de- 
, fined "humps". The first one - peaking somewhere between the IR 

and the X-ray bands - derives from the synchrotron emission of 
relativistic electrons (or pairs) in the jet. The origin of the second 
component which exhibits a maximum at 7-ray energies is more 
debated. Leptonic models attribute it to the inverse Compton emis- 
sion of the same electrons responsible for the lower energy bump. 
Hadronic models, instead, assume that the 7 rays are the leftover 
of reactions involving relativistic hadrons. Blazars are further di- 
vided into two broad groups, BL Lac objects and flat spectrum ra- 
dio quasars (FSRQs). BL Lacs are defined by the weakness (or even 
absence) of thermal features (most notably broad emission lines) in 
their optical spectra. This evidence leads to the common belief that 
the nuclear region of BL Lacs, where the jet forms and accelerates 
is rather poor of soft photons. On the contrary, FSRQs display lu- 
minous broad (t; > 1000 km s"^) emission lines, indicating the ex- 
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istence of photo-ionized clouds rapidly rotating around the central 
black hole and organized in the so-called broad line region (BLR). 

Besides their importance for to the study of the structure and 
functioning of relativistic jets, growing interest for blazars is moti- 
vated by the use of their intense 7-ray beam as a probe of the extra- 
galactic background light (EBL) (see e.g. Aharonian et al. 2006) 
and the intergalactic magnetic fields (see e.g. Neronov & Vovk 
2010, Tavecchio et al. 2010) and - even more fundamentally - for 
the study of new physical phenomena beyond the Standard Model 
(SM) of elementary particle interactions and even quite radical de- 
partures from conventional physics such as violation of Lorentz in- 
variance (for a review, see Liberati & Maccione 2009). Concerning 
the possible evidence for new elementary particles, particular im- 
portance is played by axion-like particles (ALPs) - namely very 
light spin-zero bosons a characterized by a two-photon coupling 
a77 - which are a generic prediction of many extensions of the 
SM (for a review, see Jaekel & Ringwald 2010). In the presence 
of an external magnetic field, the 077 coupling gives rise to the 
phenomenon of photon-ALP oscillations (conversion), quite simi- 
lar to oscillations of massive neutrinos of different flavour. Within 
the present context, photon-ALP oscillations can drastically alter 
the propagation of hard photons of energy E from a blazar to us. 
Two different scenarios have been proposed. One - called DARMA 
- contemplates photon-ALP oscillations 7 — a — s- 7 as taking 
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place in intergalactic space (De Angelis et al. 2007, De Angelis 
et al. 2009, Mirizzi & Montanino 2009, De Angelis et al. 2011 
[D2011]) where a large-scale magnetic field in the nG range is 
supposed to exist, which is consistent with present upper bounds 
(Kronberg 1994, Blasi et al. 1999, Grasso & Rubinstein 2001) and 
AUGER results (Abreu et al. 2010). Alternatively, the conversion 
7 — >^ a can occur inside the blazar jet whereas the re-conversion 
a — !> 7 can happen in the Milky-Way (Simet et al. 2008). Needless 
to say, also both options can be realized at once (Sanchez-Conde 
et al. 2009). In either case, what we are used to simply regard as 
a photon behaves for some time as a "true photon" and for some 
time as an ALR Now, "true photons" can disappear from the beam 
along their way to us through the 77 — >■ e^e" scattering with low- 
energy photons of the EBL while ALPs propagate unaffected by 
the EBL, since the process 07 ^ 7 is kinematically forbidden. As 
a result, owing to photon- ALP oscillations the effective photon op- 
tical depth Tcs{E) gets smaller, and even a tiny decrease entails a 
large enhancement of the photon survival probability P-,^-, (E) as 
compared with conventional physics since 

Pj^j{E) = exp[-reff(£)] . (1) 

In addition, because the EBL dimming increases with E whereas 
the photon-ALP oscillations alone are i5-independent, the result- 
ing observed blazar spectra become harder than expected and in 
particular the "7-ray horizon" gets considerably enlarged. 

Coming back to blazar observations, the recent evidence of 
VHE emission by FSRQs poses a quite serious challenge. As 
stressed above, the surrounding of the inner jet in FSRQs is rich 
of optical-UV photons emitted by the BLR, necessarily implying 
a huge optical depth for 7 rays above 10-20 GeV (see e.g. Liu & 
Bai 2002, Tavecchio & Mazin 2009 [TM2009], Poutanen & Stem 
2011). Therefore, the observation of some FSRQs at TeV energies 
raised great surprise (Albert et al. 2008, Aleksic et al. 201 la, Wag- 
ner & Behera 2010). Moreover, the detection of an intense VHE 
emission from PKS 1222-1-216 (Aleksic et al. 2011a) observed by 
MAGIC to double its flux in only about 10 minutes - thereby flag- 
ging the extreme compactness of the emitting region - is very 
difficult to fit within the standard model (Tavecchio et al. 2011 
[T2011]). 

So far, the only possibility to solve the apparent contradiction 
between the detection of the intense and relatively hard flux and 
the rapid variation appears to invoke the existence of very small 
{r ^ 10^^ cm) emitting regions beyond the BLR {R ~ 10^** cm), 
that is at a large distance from the central engine. 

Our proposal to explain the observations of PKS 1222-1-216 as- 
sumes the validity of the standard model for photon production but 
we add the new assumption that photon-ALP oscillations take place 
inside the source in such a way that a considerable fraction of pho- 
tons can leave it, indeed much in the same fashion that a consider- 
able fraction of photons emitted by blazars avoid EBL absorption in 
intergalactic space. We envisage that the 7 — > a conversion occurs 
before most of the photons reach the BLR. Accordingly, ALPs tra- 
verse this region unimpeded and outside it the re-conversion 0^-7 
takes place either in the same magnetic field of the source or in that 
of the host galaxy. So, our proposal differs from any previously 
considered one. 

The plan of the paper is as follows. The mechanism of photon- 
ALP oscillations is reviewed in general terms in §2. Our proposed 
model for PKS 1222-1-216 is described in §3 and its results are pre- 
sented in §4. Finally, in §5 we compare our model with the two-blob 
scenario developed in [T201 1] to explain the same observations and 
we discuss some features of the model presented in this Letter. 



2 PHOTON-ALP OSCILLATIONS 

We consider in the following a monochromatic photon beam of en- 
ergy E travelling along the y direction from the centre of the source 
to us. In this section we use natural units, in which c = h — 1. 

Below, we cursorily describe the mechanism whereby photon- 
ALP oscillations occur in the beam propagating in a generic ionized 
and magnetized medium where photons can be absorbed but ALPs 
can not (a very detailed account can be found in 

[D201lf])- 

Basically, the SM is enlarged so as to include an ALP which 
is described by the Lagrangian 

-Calp = i d'^ad^a - i + E • B a , (2) 

where m is the ALP mass and M is a constant with the dimension 
of an energy. More specifically, it is generally supposed that m < 
1 eV, and the CAST experiment (Arik et al. 2009) sets the robust 
bound M > 1.1 X 10"' GeV for m < 0.02 eV while from the 
absence of 7 rays coming from SN 1987 A the stronger bound M > 
10^^ GeV has been derived for m < 10"^" eV even if with a large 
uncertainty (Brockway et al. 1996, Grifols et al. 1996). Clearly the 
last term in Eq. ^ gives rise to the ajj coupling mentioned above. 
Note that because of the specific form of this coupling oc E ■ B /M 
only the component of B in the x — z plane By couples to a, 
and in addition all physical results depend only on the combination 
Bt/M. 

The first point that deserves concern is that in the approxima- 
tion E ^ m the beam propagation equation becomes (Raffelt & 
Stodolski 1988) 

\^i — +E + M{y)j I A,{y) 1=0, (3) 

where Ax{y) and ^2(2/) are the two photon linear polarization 
amplitudes along the x and z axis, respectively, a{y) denotes the 
ALP amplitude and M represents the photon-ALP mixing matrix, 
whose explicit form is (Mirizzi et al. 2005) 

/ A,, AS^ \ 
M{y) =\ A22 Al^ . (4) 
\ AS^ Al., A,„ / 

Denoting by u)pi{y) the plasma frequency - related to the 
electron number density ne{y) by (jJp\{y) — 3.69 x 
10^^^ \/ ne{y) /cm'-' eV - by X^{E, y) the photon mean free path 
and by Bx{y) and Bz{y) the x and z components of the magnetic 
field B(j/), the various A terms are Aaii; = A^z = —L0pi{y)/2E + 
i/2X^{E,y),Al^ = B4y)/2M,Al^ = B,{y)/2M mdAaa = 
—m^ /2E. A very remarkable fact is that Eq. l|3j is a Schrodinger- 
like equation with y playing the role of time, and so the relativis- 
tic beam in question can formally be treated as a three-level non- 
relativistic unstable quantum system with non self-adjoint Hamil- 
tonian H{y) — ~[E + Ai{y)]. This circumstance allows us to 
describe the behaviour of the beam by non-relativistic quantum 
mechanics. Moreover, we denote by U{y, yo) the transfer matrix, 
namely the solution of Eq. ^ with initial condition U {yo, yo) = 1- 
Actually, we assume that initially the beam is fully made of pho- 
tons whose polarization is unknown to a considerably extent, so 
that it is safe to regard the beam initially as unpolarized. There- 
fore, it must be described by the (generalized) density matrix p{y) 

^ Although the results presented in this paper are all correct, Eqs. (104) 
and (134) are written incorrectly and their denominator should be dropped, 
so that they have to look like Eq. (48). 
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which - thanks to the above analogy - obeys the Von Neumann-hke 
equation 

= pM'' -Mp (5) 

dy 

associated with Eq. Then it follows that even ifM^'T^M the 
solution of Eq. ^ is given by p{y) = U{y,yo) p{yo) U^{y,yo), 
and since in the 7-ray band the photon polarization cannot be mea- 
sured the photon survival probability after a distance y for the beam 
initially in the state unpolarized state punp is 

P-,^-,{E,y)= J2 Tr(p.t/(y,0)pu,pf/t(y,0)) , (6) 

i — x^z 

where p^ and pz denote the linearly polarized states along x and z, 
respectively. 

Consider now the particular case in which no absorption is 
present and B is homogeneous. Then we can choose the z axis 
along Bt so that Bx ~ 0. In this case the photon-ALP conversion 
probability can be computed exactly and reads 

where we have set A^^^ = [(m^ - ujli)/{2E)f + [Br/Mf. 
Hence, we see that for E Et with E, = \m? -u}'^^\M/2Bt the 
strong-mixing regime takes place, in which the photon-ALP con- 
version probability becomes maximal and energy-independent. An- 
other important consequence of the strong-mixing regime is that the 
general mixing matrix A^(y) simplifies since the terms i^pi(j/) /2E 
and m? /2E should be discarded. 



3 A MODEL FOR PKS 1222+216 

We now proceed to build up our model for PKS 1222+216. 

In the first place we need the knowledge of the relevant phys- 
ical parameters. We assume a disk luminosity Ld — 5 x 10*^ erg 
s~\aradius of the BLR iJsLR 7xl0" cm ([T20 11]), and stan- 
dard values of the BLR cloud number density ric — 10^" cm"'' and 
temperature Tc ~ lO" K (see e.g. [TM2009]). Since the filling fac- 
tor of the clouds is small, the average electron number density rie 
relevant for the beam propagation is much smaller than Uc. Mod- 
els attributing the confinement of the clouds to a hot — 10^"* 
K external medium in pressure equilibrium with the clouds yield 
He ~ lO^"'^ cm~^ (see e.g. Krolik et al. 1981). However, the pres- 
ence of such a hot confining medium is disfavored by the lack of 
the necessarily expected bright X-ray emission. So, an extra con- 
tribution to the pressure confining the BLR clouds is expected e.g. 
from the magnetic field (see e.g. Netzer 2008). What matters for us 
is that Tie gets considerably reduced, and for definiteness we take 
ne ~ 10^ cm"''. At any rate, since we are working in the strong- 
mixing regime the photon survival probability is unaffected by rie, 
which only enters £*. What about the magnetic field B? Unfortu- 
nately, its strength and geometry in the innermost region of an AGN 
are very poorly known and in first approximation B can be taken 
as uniform. A possible estimate of its strength can be derived by 
assuming equipartition between magnetic and gas pressure, lead- 
ing to i? ~ 1 G (see e.g. Rees 1987). Thus, we assume B ~ 1 G 
from the emission region out to -Rblr. But since we need to follow 
the beam until it exits from the host galaxy we have also to provide 
a description of B beyond the BLR. Due to the fact that the BLR 
is surrounded by ionized gas with a smoothly declining density, it 
looks natural to suppose that the B field lines are frozen-in, so that 



flux conservation entails B{R) ~ (-R/-Rblr)~^ G. Furthermore, 
FSRQs are hosted by elliptical galaxies, whose B is very poorly 
known. Nevertheless, it has been argued that supernova explosions 
and stellar motion give rise to a turbulent B ranging from about 
10 pG at the centre to about 1 pG at 10 kpc (Moss & Shukurov 
1996). Because of such a small gradient of B across the galaxy and 
its associated uncertainty, we will take B ~ 5 pG throughout the 
galaxy. 

We start by computing the optical depth t{E) of the beam 
photons according to conventional physics, which amounts to con- 
sider only the process 77 — >■ e'^e^ . We follow the same procedure 
developed in [TM2009]), to which we refer the reader for a full 
description of the method. The optical depth is derived using the 
detailed spectrum emitted by the photo-ionized BLR clouds. The 
result is plotted as the blue dashed line in Fig.[T] 

Next, we turn our attention to the situation in which photon- 
ALP oscillations occur in the beam. Owing to the varying physical 
properties in the space traversed by the beam, we perform an inde- 
pendent analysis of three regions: 

• The inner part of the blazar, namely from the centre to 7?blr, 
where B is supposed to be uniform. 

• The region surrounding the BLR, where B is assumed to de- 
crease like R^^ until it reaches the strength of the galactic B at the 
galactocentric value _R* ~ 10^ pc. 

• The galactic region from _R* ~ 10^ pc to iihost — 10 kpc 
where B has a constant strength but a turbulent nature. 

We begin from the first item. The most important issue to settle 
is to find out the actual value of i5, above which the strong-mixing 
regime sets in for the assumed values of the electron density and of 
the strength of B. To this end, it is useful to rewrite E, in the form 

i...0.26xl0-|^^-.^,|(^) (^^) eV-.(8) 

For Tie ~ 10'' cm"'' we get ujpi ~ 3.69 x 10"'"* eV. Moreover, we 
have chosen B ~ 1 G but we do not know its direction, which is 
quite likely to change with the distance from the centre. Because 
the photon-ALP oscillation vanishes if B is along the beam while 
it is maximal for B normal to the beam, we suppose that B is on 
average at an angle of 45° with the beam direction. Therefore we 
have Bt ~ 0.7 G. Finally, the bound M > 1.1 x 10^° GeV cer- 
tainly holds. So, putting everything together, from Eq. ([Sj we get 

~ 0.2 MeV under the assumption that m < ujpi [D2011]. 
This is a welcome conclusion, because we do not run the risk that 
photon-ALP oscillations affect X-ray observations of FSRQs while 
they do affect their observations performed both by Fermi/LAT and 
by Cherenkov telescopes like MAGIC. 

Moreover, we get \-y{E) — RB-Li\./r{E) for R < i?BLR. 
At this point, following [D2011] the calculation of the transfer 
matrix - which we denote by ?7i(-Rblr,0) - is in the present 
case straightforward. As a consequence, the beam state at -Rblr 

is PbLR = (7l(i?BLR,0)punp(7i(-RBLR,0). 

Let us next turn our attention to the second item. Remarkably 
enough, for B(i?) ~ (i?/i?BLR)"^ G Eq. ID can be solved exactly 
- in this region absorption can safely be neglected - and we denote 
by (72 (i?*, -Rblr) the corresponding transfer matrix. So, the beam 
state at Rt is pR, = (72(-R,, -Rblr) Pblr [/J(-R., -Rblr) = 
U2{R„Rbi.r) (7i(-Rblr, 0) punp UI{Rblk, 0) 17+ (-R., -Rblr). 

Finally, we address the third item. According to Moss & 
Shukurov (1996), we model B in the host elliptical galaxy as a 
domain-like structure, with strength 5 pG and domain size equal to 
150 pc. In this case the beam propagation can be computed ex- 
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4 RESULTS 




10 100 
E [GeV] 

Figure 1. 77 optical depth as a function of energy for photons propagating 
in the BLR of PKS 1222+216 (blue dashed line), calculated within conven- 
tional physics as in [T2009] . The red solid line shows the effective optical 
depth assuming photon- ALP conversion. See text for details. 
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Figure 2. Photon survival probability Pj^j at 1 GeV (blue dashed line) 
and at 200 GeV (red solid line) as a function of the inverse a-yy coupling 
constant M. The gray vertical line indicates the value used in the model, 
M = 4 X 10" GeV. 



Our source has been observed by FermifLAT in the energy range 
0.3-3 GeV and by MAGIC in the band 70 - 400 GeV. There- 
fore we focus our attention on the energies E — 1 GeV and 
E — 200 GeV as representative of the two kinds of measurements. 
We show in Fig. |2] P-y_>-, (E) at these energies as a function of 
M, namely the inverse ayy coupling constant: The blue dashed 
line corresponds to 1 GeV while the red solid line to 200 GeV. 
We see that P^_>.y(200 GeV) reaches its maximum close to M = 
4 X 10^^ GeV, which is evidently our most favourable case. Since 
at 1 GeV no photon absorption occurs, on the basis of Eq. l|7j we 
expect an oscillatory behaviour, even though outside -Rblr B first 
decreases and then has a random domain-like structure. Moreover, 
since the oscillation wavelength is Lose ~ 2tiM /B it has to in- 
crease with M. Both features are reproduced by the blue dashed 
line. The behaviour of the red solid line is more difficult to under- 
stand in intuitive terms due to the presence of a strong photon ab- 
sorption. Certainly beyond a certain value of M photon- ALP oscil- 
lations become ineffective and absorption dominates: we see from 
the figure that this is indeed the case for AI > lO^'' GeV. For 
smaller values of M photon- ALP oscillations become efficient, but 
when they are very efficient - and this occurs fairly soon ~ some 
ALPs are re-converted back into photons before leaving the BLR, 
thereby getting absorbed. This explains why P^^^{E) peaks at 
some value of AI. 

In order to get deeper insight into the physical relevance of our 
result, using Eq. ([T} from P^^^ {E) we obtain the effective optical 
depth Tcff (i?) - as evaluated for Af = 4 x lO" GeV - which is 
represented by the red solid line in Fig.[T] The effect of the photon- 
ALP oscillations on the beam propagation can be readily appreci- 
ated. Indeed, photon- ALP oscillations lead to a drastic reduction of 
the optical depth in the optically thick range: in the MAGIC band 
it is almost constant Toff ~ 4, corresponding to a survival probabil- 
ity of about 2% as it is evident also from Fig. |2] On the contrary, 
in the optically thin region below ~30 GeV, the optical depth in 
the presence of the photon-ALP conversion is larger than the stan- 
dard one, which instead goes rapidly to zero below 10 GeV. This 
behaviour can be simply understood: a fraction around 40% of the 
7 rays originally emitted by the source and converted into ALPs 
do not reconvert to photons, therefore leading to a reduction of the 
observed photon flux. 

Let us now explicitly address the impact of our model for 
the SED of PKS 1222-F216, which is illustrated in Fig. |3] The 
red points at high energies are the FermifLAT and MAGIC EBL- 
deabsorbed data around the epoch of the VHE detection taken 
from Aleksic et al. 201 la (see [T201 1] for details and references). 
As far as our model is concerned, the intrinsic flux Fint{E) is 
clearly related to the EBL-deabsorbed observed one FohsiE) by 
Fint{E) = Fob,{E)/P^^^{E) = Fobs{E) exp[Teff (£;)]. Appli- 
cation of this relation to the red points yields the corresponding 
black points in Fig. |3] We see that the black points corresponding 
to FermifLAT and MAGIC observation describe a Compton bump 
peaking at ~ 50 GeV, similar to that of BL Lacs. 



actly as in [D2011] without absorption, and the resulting trans- 
fer matrix is denoted by U3{RhoBt, R*)- However, this last ef- 
fect turns out to give a very small contribution. Thus, by replac- 
ing W(- • ■) ^ U3{Rho.t,R,)U2{R„RBLR)Ul{RBLR,0) in Eq. 

(134) of [D2011] and going through the same steps we find the 
photon survival probability P^-^~f{E) from PKS 1222-1-216 to us. 
We stress that this conclusion ultimately rests upon the application 
of Eq. 



5 DISCUSSION AND CONCLUSIONS 

We find it instructive to compare our proposal with the one pre- 
sented in [T201 1], where the SED comprising the red 7-ray data is 
reproduced by a model consisting in two active compact regions, 
one responsible for the emission from IR to X-rays and the other, 
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Figure 3. SED of PKS 1222+216 around the epoch of the MAGIC detec- 
tion. Red points at high and very high energy are the spectrum recorded by 
FermifLKI and MAGIC (EBL-corrected). The black points represent the 
same data once corrected for the photon-ALP oscillation effect. The curves 
show the result of a standard leptonic emission model assuming two re- 
gions located inside the BLR responsible for the emission at low and high 
energy, respectively. The gray data points below 10^" Hz are irrelevant for 
the present discussion (details can be found in [T2011]). 



very compact, accounting from the rapidly varying 7-ray emission. 
The latter region is assumed to lie at J? > lO'^* cm., namely be- 
yond the BLR, so as to avoid the resulting huge photon absorp- 
tion. We refer to [T201 1] for a full discussion of the problems and 
the detailed description of the model. Briefly, each region is speci- 
fied by its size r, magnetic field B, bulk Lorentz factor F, electron 
normalization K, minimum, break and maximum Lorentz factors 
7min, 7b, 7max and slopes ni, n2. The electrons radiate through 
synchrotron and inverse Compton emission (considering both the 
internally produced synchrotron radiation and the external radia- 
tion of the BLR). For the larger region we use the parameters re- 
ported in [T201 1] while for the compact 7-ray emitting region the 
model is specified by the following parameters: r = 2.2 x 10^* 
cm, B = 0.025 G, F = 17.5, if = 4.1 x 10^ 7^1^ = 1.5 x 10^ 
7b = 10^, 7max = 4 x 10^ and slopes n\ = 2.1, n2 — 4.2. As 
discussed in [T2011] the relative position of the two regions is not 
relevant for the emission properties. 

We remark that the emission model just discussed is merely 
an example, and different and possibly more realistic scenarios can 
be constructed. The key point we want to make is that with the 
photon-ALP oscillation mechanism at work the emission can well 
originate inside the BLR just like in conventional models. 

Needless to say, our scenario naturally applies also to the other 
FSRQs detected at VHE like 3C279 and PKS 1510-089 (Albert et 
al. 2008, Aleksic et al. 201 lb, Wagner 2010), although these cases 
appear less problematic for the external emission scenario due to 
the absence of evident rapid (t < 1 day) variability. 



Finally we would like to mention a possible test of the pro- 
posed model. Because photon-ALP oscillations can mitigate - but 
not completely avoid - the 7-ray absorption inside the BLR, a nat- 
ural prediction is that at the optically-thin/optically thick transition 
around 30 GeV (in the source frame) the spectrum should display 
a feature. So, the absence of such a feature would be hard to ex- 
plain in our model but would directly support scenarios in which 
the emission occurs outside the BLR. 

We find it quite tantalizing that precisely the most favourable 
value A/ ~ 4 X 10^^ GeV for the effect considered in this Letter 
coincides with the most favourable value in the DARMA scenario 
that enlarges the "7-ray horizon" and provides a natural solution to 
the cosmic opacity problem [D2011] (also a very small ALP mass 
m < 10^" — lO"^'' eV is common to both models). Yet, we have 
neglected photon-ALP oscillations in intergalactic space - which 
would increase P^-*^ (E) ~ in order to put ourselves in the less op- 
timistic case. Remarkably enough, an independent laboratory check 
of our proposal can be performed by the planned upgrade of the 
ALPS experiment at DESY. 

A much more detailed presentation of the matter reported here 
will soon appear. 
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